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Possibility of a novel pseudo-scalar (octupole) order is studied theoretically for orbitally de- 
generate systems with strong spin-orbit coupling such as Ce^^Lai-i^Be. It is discussed that coex- 
istence of an octupole order parameter and antiferromagnetic fluctuation should lead to drastic 
softening of the elastic constant by a mode-mixing effect. Nonlinear coupling between dipole, 
quadrupole and octupole fluctuations is taken into account in terms of a Ginzburg-Landau-type 
functional which is derived microscopically through path integral. 
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In this paper we study possibility of a new pseudo- 
scalar order in orbitally degenerate systems with strong 
spin-orbit interaction. Our motivation is a strange phase 
called IV found in Cea;Lai_a;B6 with x ~ 0.75. The phase 
IV has the following properties: 

(i) The transverse ^elastic constant C44 shows a drastic 
softening (~ 20%). B) 

(ii) There is almost no magnetoresistance in contrast 
with other phases H' 

(iii) The phase IV is isotropic magnetically in contrast 
with phases II and III.eP 

Although the magnetic susceptibilityj-, shows a cusp at 
the phase transition to the phase IV,H^ preUminary ex- 
periment of neutron scatterin^^ has found no magnetic 
Bragg scattering along high symmetry axes such as (100), 
(110) and (111). The experimental fact (iii) suggests 
that the order parameter is a scalar or a pseudo-scalar 
instead of a dipole or a quadrupole. A pseudo-scalar or- 
der may bring a gap or a pseudo-gap in the magnetic 
excitation spectrum with zero momentum since there is 
no Goldstone mode in contrast with the Neel order. It 
is conceivable that this is related to the fact (ii) . 

These facts together with the most conspicuous fact (i) 
suggest that a new type of order is involved in the phase 
IV. We hypothesize that the pseudo-scalar component of 
the octupole moment is the order parameter of the phase 
IV, and explore some consequence of the hypothesis. In 
particular we provide possible mechanism to explain the 
fact (i) in terms of the octupole order. 

In order to characterize the new order we begin with 
the symmetry analysis of the Fg level which is the 
crystalline-electric-field ground state of each Ce^+ ion 
in Ggx^&i-xQq- The excited level Fy lies about 500 K 
from the Fg level, and can be safely neglected for our 
purpose. The four states in the Fg level are represented 
with use of the basis | J^) of J = 5/2 states as follows: 



The four-fold degeneracy can be broken in a variety of 
ways depending on types of symmetry-breaking fields: 

(a) magnetic order — In this case there remains no de- 
generacy since the Zeeman splitting caused by internal 
magnetic field is different between the two orbitals. 

(b) orbital order — The two-fold Kramers degeneracy 
remains since the time-reversal symmetry is not broken. 
In the phase II of CeBg the antiferro-orbital order (or 
the quadrupole order) is realized. In the phase III simul- 
taneous presence of (a) and (b) have been invoked for 
explanation of the neutron scattering experiment. One 
can then ask whether it is possible to have 

(c) breakdown of the time-reversal symmetry without 
breaking the orbital degeneracy. 

In this paper we assert that the question can be answered 
in the affirmative. 

To describe the spin and orbital degrees of freedom 
in a concise way, we introduce the pseudo-spin opera- 
tors which are represented by two sets of Pauli matrices: 
{gx, cTy^az} and {tx, Ty,Tz}. The former operates on the 
Kramers partners, and the latter on the orbital partners. 
Namely we have 

with a = 1,2. The x and y components of pseudo-spins 
change from one state to another in the Fg subspace. 
We can express the physical operators adapted to the 
point-group symmetry pU^ng the pseudo-spins. They are 
enumerated as follows:ErLP 
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where we have introduced Hnear combinations of and 
as 

ry± = i(±v^r--r^), (3) 

C± = -i(r-±V3r^). (4) 

The subscript u represents the odd property under the 
time reversal, and g does the even one. 

Among these operators, in V-^u has the same matrix 
element as the symmetrized product of J^^J^J^', and is 
regarded as a component of the octupole moment ten- 
sorou-* This component commutes with discrete rota- 
tions of the cubic group, and is odd under the time re- 
versal. The latter property is consistent with being 
pure imaginary. Hence is regarded as a pseudo-scalar 
operator. On the other hand the operators belonging to 
T^g and T^g describe the quadrupole operators which are 
even under the time reversal. The dipole operators un- 
der the cubic symmetry is decomposed into r4^^ and v'f^ . 
For example the x-component of the magnetic moment 
is given in units of the Bohr magneton by 

Af- = a- + ^7?+a^ (5) 

where the first term on the right-hand side belongs to 
r4^^ and the second one to T^l^ . We note that each three- 
dimensional odd representation is a linear combination 
of dipole and octupole operators. In other words, dipole 
and a part of octupole operators mix under the point- 
group symmetry. The remaining representation T^^ cor- 
responds to pure octupole operators other than , and 
describes a part of third-rank tensors composed of J" . 

We take the simplest possible model to describe the 
coupling among dipole, quadrupole and octupole mo- 
ments. The conduction electrons which give rise to the 
Kondo effect are not treated explicitly. For the magnetic 
moment, we keep only the representation V^}^ and ne- 

(2) 

gleet the orbital dependent part for simplicity. The 
effect of the neglected part will be discussed in the end of 
the paper. For the quadrupole moment, we include only 
the Fsg component which is known to be dominant in 
CeBg. In addition to these we include the pseudo-scalar 
component of the octupole moment. The Hamilto- 
nian is given by 

i¥=3 

M . .... 

where j\j with a = m^e,8 describes intersite interac- 
tion of either magnetic, electric (quadrupole) or octupole 
degrees of freedom. 

Since the quadrupole operator r^cr is the product of 
the dipole cr and octupole operators, there is intrinsic 
frustration if all of them favor the antiparallel arrange- 
ment between nearest-neighbor pairs. The complicated 
magnetic structure of phase III seems to be realized as 
a compromise between these competing interactions. If 
there is a delicate balance in realizing the actual struc- 
ture, slight change of the balance with substitution of 
Ce by La might lead to another structure. Although a 



quadrupole order occurs first as temperature is lowered 
in pure CeBg , substitution of La may favor a breakdown 
of the time-reversal as the first instability. We interpret 
the phase IV as being realized in this way. 

The pseudo-spin representation is introduced just to 
reproduce the matrix elements of multipole operators. 
Hence any approximation to decouple a and t is not 
meaningful physically. Instead we should consider on 
equal footing each fluctuation with a point-group sym- 
metry. To deal with coupled multipole fluctuations 
from the high-temperature side, we work with the path- 
integral representation of the partition function. Namely 
we use the Stratonovich-Hubbard identity to replace the 
intersite interaction by the local interaction between aux- 
iliary fields and multipole moments. The most impor- 
tant technical point is that we can translate the full non- 
linearity represented by the pseudo-spin operators into 
the explicit coupling term in the Ginzburg-Landau-type 
(GL-typc) functional. After this translation we can in- 
troduce a suitable approximation such as the mean-field 
theory with respect to these auxiliary fields. 

The local interaction at each site i is of the form 

(t^^-a^+ir^^rf + ^^. (7) 

These auxiliary fields 4>^ , cr j and obey the Gaussian 
distribution. In this paper we confine ourselves to the 
static approximation where dynamical fluctuation is ne- 
glected. Although some quantum effects already escape 
at this stage, our approximation keeps faithfully the non- 
trivial commutation property of pseudo-spins. Hence 
we expect that interesting consequences due to coupling 
between different multipoles can be understood qualita- 
tively within the static approximation. 

In computing the the partition function, we first carry 
out the trace over orbital part taking such basis that 
makes r** diagonal. Then we are left with another trace 
over the Kramers partners. We note that an effective 
magnetic field which couples with cr^ is given by (p^ ± 
^j. Then the latter trace is performed most easily by 
rotating tkp-jquantization axis so that ± is along 
the z-axiga-llJ) of cri. In this way we obtain 

Z = Trexp(-/3i7) = j VcpViViP cyi]){~f3T), (8) 

where the functional T consists of three parts: T = 
—NThiA + Tq + J- 1. The first part is the entropy term 
and the second one, which describes the Gaussian fiuc- 
tuation, is given by 

^ q 

[Jeiq)-' - mq ■ ^ q + [Jsiq)-' - /3]V^qV^-q}, (9) 

where Ja{q) with a — to, e, 8 is the Fourier transform 
of Jif' ■ On the other hand the interaction part Ti is 
derived in the closed form from 

exp[-/3(.Fi + T,n)] = [] 2[cosh(/3|c/,, + exp(-/3V») 

i 

+ cosh(/3|0, -|J)exp(/?^,)], (10) 
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where !Fq2 consists of the entropy term and the second- 
order part which is already accounted for as the terms 
proportional to /? in eq.(9). By expanding !Fi up to 
fourth order with respect to auxiliary fields we obtain 



12 



This GL-type functional provides us with the starting 
point for discussing growth of the order parameter and 
coupling among magnetic, quadrupole and octupole fluc- 
tuations. 

The relation between the octupolar susceptibility 
Xsig) and the fluctuation (i^qijj-q) is given by 



(12) 



We have analogous relations also for the magnetic sus- 
ceptibility Xm{q) and for the quadrupolar susceptibility 
Xe{q)- If one neglects !Fi in taking the thermal average, 
one obtains the RPA result given by 



Xa{q) = f3/[l- Js{q) 



(13) 



Let us assume that the high-temperature phase becomes 
unstable first against formation of an octupole order with 
the wave vector Q. The transition temperature Tg, is de- 
termined in the RPA (or the mean-field approximation) 
as 



Ts = JaiQ). 



(14) 



One could include fluctuation corrections coming from 
J^i . However such sophistication is not necessary to our 
purpose of demonstrating the mode mixing in the pres- 
ence of an octupole order. 

The magnitude of the order parameter (V'q) can be 
determined by the mean-field theory which reduces to 
the standard GL theory near the transition temperature. 
We assume that Q is given by (1/4, 1/4, 1/4) in units of 
27r/a where a is the lattice parameter of the cubic crystal. 
This choice of Q is motivated by the known magnetic and 
quadrupole patterns in the phase IIlJj' Namely the mag- 
netic supercell in the (001) plane contains eight Ce ions 
in the 2^/2 x 2\/2 structure and orbital supercell contains 
two Ce ions in the \/2 x V2 structure. To be consistent 
with these structures the octupole order should also have 
the 2\/2 X 2\/2 structure which corresponds to the wave 
vector (1/4, 1/4) in the plane. Assuming the cubic sym- 
metry at Tg we expect the wave vector (1/4,1/4,1/4) 
(together with its stars) as a reasonable candidate. 

In the mean-field theory we replace fluctuating fields 
in T by their averages. Then the stationary condition 
5T /6{'iljQ) — leads to the result 



(15) 



where N is the total number of unit cells, and we take 
into account the fact that there are eight equivalent Q's. 
In ultrasonic measurement, the external strain is cou- 



pled to quadrupole moments by 

Hext = ^^g(r)er7(-R*)Or7(-R0, (16) 



r7 



where g(r) is the coupling constant of a representation 
r, which runs over Tzg and Fsg, and the strain er^yiRi) 
at site i and the quadrupole moment Orj{Ri) belong 
to the same representation. For the quadrupole moment 
with the T^g symmetry, the components Or-y correspond 
to Cartesian ones Oa/s as 



T%^,r^a^) 



(17) 



The change AC44 of the elastic constant and the change 
Avs of the transverse sound velocity are related to the 
homogeneous quadrupolar susceptibility Xe by 

5(r5,)^ 



ACa 



2Av, 



-Xe 



(18) 



C44 Vs Mvl 

where M denotes the mass of the unit cell. Therefore 
the enhanced Xe leads to softening of the C44 mode. 

The T^u order breaks the time-reversal invariance but 
leaves the orbital degeneracy in contrast with the mag- 
netic order. With a finite octupole order parameter, 
the orbital fluctuation hybridizes with the magnetic fluc- 
tuation. Then the quadrupolar susceptibility measured 
by ultrasound probes the magnetic fluctuation with the 
wave vector Q. Note that the wave number of ultra- 
sound is negligible as compared with |Q|. We consider 
the z-component of </)q and ^q, and omit writing the 
component index since other components follow the same 
equation. In the RPA, we obtain the following equations 
to determine Xe'- 



f3 



Je{OfXe^P{il)~Je{0), 



(19) 



(20) 



where we take the order parameter real. This equation 
applies to the temperature range where neither magnetic 
nor quadrupole order is present. Solving eq.(20) we ob- 
tain 



Xe 



1 



-Jr. 



Je{0) [{1 - Jm)il ~ Je) ~ JraJeP^^JQ)' 



1 



where j,n ~ Jm[Q)P and je = Je(0)/3. It is easily seen 
that the result reduces to the conventional mean-field 
one above the transition temperature Tg. 

Figure 1 shows an example of numerical results with 
tentative values of interactions: Jm{Q) /Tg, — 0.6 and 
Je(0)/r8 = —0.2. We have taken the negative value 
of Je(0) to be consistentent with the nearest-neighbor 
antiferro-quadrupolar interaction. It is seen that the 
quadrupole susceptibility increases significantly below 
Tg. The reason for the increase is the increased coupling 
with growing antiferromagnetic correlation with the 
wave vector Q. Note that this antiferromagnetic correla- 
tion is closely related to the wave vector (1/4, ±1/4, 1/2) 
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of the magnetic order in the phase III. Hence the growth 
of this correlation is naturaUy expected. 
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Fig. 1. The quadrupolar susceptibility plotted as a function of 
temperature. The energy is normalized by Tg = Js{Q)- The 
interactions are taken to be Jm{Q) = 0.6Tg and Je{0) = — 0.2T8. 



With the wave vector Q = (1/4, 1/4, 1/4) the transi- 
tion between the phases IV and III should be of first 
order since the phase III has different wave vectors 
(1/4, ±1/4, 1/2) for the superstructure. This is consis- 
tent with experimental observation^) The transition to 
the phase II with increasing magnetic field should be also 
of first order with the present Q, since the phase II has 
the wave vector (1/2,1/2,1/2) as the quadrupole order, 
which should mix with the octupole order in magnetic 
field. Experimentally the phase boundary is extremely 
narrow; the boundaries to phases I and III seem to merge 
with the II-IV boundary in the (T, H) plane with H par- 
allel to (001). Our model is probably too simple to be 
applied to such details of the phase diagram. If Q of the 
octupole order were the same as that of the quadrupole 
order, the phase boundary Jbetween the phases II and 
IV should have disappearedO This is similar to the ab- 
sence of phase boundary between the Neel phase and the 
antiferro-quajdrupolar phase in TmTe under finite mag- 
netic fieldMe 

Recently a phase diagcap analogous to that of 
Cea;Lai_a;B6 has been foundlij' in a tetragonal compound 
II0B2C2 for magnetic field along (110). It appears that 
the phase IV in this compound has a Neel order. It 
should be interesting to see whether there is an elas- 
tic anomaly in the phase IV of II0B2C2. We note that 
the angular momentum of the Hund-rule ground state 
of Ho'^+ is as large as J = 8, and that 4f electrons here 
are more localized than in Ce^^Lai-xBg. Hence quantum 
fluctuations should be less significant in H0B2C2. 

In this paper we have proposed the simplest theory 
that can describe the coupling effect between the oc- 
tupole order and dipole as well as quadrupole fluctua- 



tions. The drastic softening of C44 is interpreted as a 
consequence of the coupling effect. We now discuss pos- 
sible directions of further development. One can com- 
pute the magnetic susceptibility Xm in a manner similar 
to what we have done for Xe- The homogeneous mag- 
netic susceptibility is influenced by the quadrupole fluc- 
tuation with the wave vector Q. It turns out that the 
susceptibility has a cusp at Tg provided that Je{Q) is 
negative. This ^behavior is similar to the one observed 
experimentally .u-* If Je{Q) is positive, on the contrary, 
Xm increases below Tg like Xe- Unfortunately we do not 
have further information on Je{Q) or Xe{Q)- The lack of 
information is in contrast with Xm{Q) which is related 
to the magnetic order in the phase HI as we discussed 
above. 

Experimentally, the ordered magnetic moment in the 
phase HI lies in the (001) plane. This anisotropy is a 
consequence of the orbital order with the Fsg symmetry 
where the wave function extends^ toward (110) or (—110) 
depending on the sublattice.BEj-' It should be possible 
to identify the octupole order if one can observe induced 
magnetic moment under uniaxial stress. Namely if the 
stress is applied along the (110) direction, we expect 
that a magnetic moment with Q is induced along the 
(001) direction. Recently large change of the magnetic 
anisotropy was found ia phases III and IV by applica- 
tion of uniaxial stress £3) Such anisotropy is taken into 
account only when one includes the component in 

addition to cr included in this paper. We plan to include 
(2) 

the component and the quantum fluctuation effects 
in a future publication. 
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